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Parameter Relations for the Shinnar-Le Roux 
Selective Excitation Pulse Design Algorithm 

John Pauly, Member, IEEE. Patrick Le Roux, Dwight Nishimura, Member, IEEE, and 
Albert Macovski, Fellow, IEEE 



Abstract— Selective excitation pulses are used in magnetic resonance 
Imaging (MRI) to isolate a specific slice through the subject. The pulse 
design problem is in general nonlinear due to the nonlinearity of the 
Bloch equation. Recently a direct solution has been proposed by Shin- 
nar and his co-workers [lJ-{6], and Le Roux [7]-[9]. 

In this paper we first present an overview of the Shinnar-Le Roux 
(SLR) algorithm. We then show how the performance of SLR pulses 
can be very accurately specified analytically. This is useful because it 
tells how to design a pulse that produces a specified slice profile. More 
importantly, it allows the pulse designer to trade-off analytically the 
parameters describing the pulse performance. We present several ex- 
amples to illustrate the more important trade-offs. These include lin- 
ear-phase and minimum- and maximum-phase pulses. Linear-phase 
pulses can be refocused with a gradient reversal, and can be used as 
spin-echo pulses. Minimum- and maximum-phase pulses have better 
slice profiles, but cannot be completely refocused. 



I. Introduction 

MAGNETIC resonance imaging (MRI) is inherently a vol- 
ume imaging modality. Spatially selective RF pulses re- 
strict this volume to a thin slab of material, which can then be 
resolved in the two remaining directions to produce a two-di- 
mensional image. The RF pulse design problem is in general 
nonlinear due to the nonlinearity of the Bloch equation. For 
small-tip-angle pulses a Fourier analysis can be used f 10], [1 1 1, 
These pulses continue to perform reasonably well for tip angles 
on the order of tt/2. Better tt/2 pulses can be designed using 
optimization methods. For tip angles of 7r, optimized pulses are 
essential for good slice profiles [12]-[16]. 

Recently a new method for slice selective pulse design has 
been proposed independently by Shinnar and his co-workers 1 1 j- 
[6J, and Le Roux [7]-[9]. It is based on a discrete approxima- 
tion to the spin domain version of the Bloch equation. The RF 
pulse design problem can then be mapped into a digital filter 
design problem, which can be solved using well-known algo- 
rithms. The resulting digital filter is then mapped back into an 
RF pulse. 

In the first part of the paper we present a brief overview of 
the Shinnar-Le Roux (SLR) algorithm. Enough details pre- 
sented to allow the reader to implement the algorithm. 

The main part of the paper concerns the specification of the 
slice profile produced by an SLR pulse. Simple analytic expres- 
sions relate the transition width, in-slice ripple, out-of-slice rip- 
ple, and time-bandwidth product. These relationships also show 
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how to design a pulse that achieves the predicted theoretical 
performance. Previously pulse sequences had to be designed 
around available RF pulses. Now RF pulse design can be an 
integral part of pulse sequence design. 

To illustrate the nature of the slice profile trade-offs we pre- 
sent several pulse design examples. The case of a slice-selective 
tt/2 pulse is considered in some detail. Then examples illus- 
trating the design of inversion, spin-echo, and saturation pulses 
are presented. 

These examples include both linear-phase, and minimum- and 
maximum-phase pulses. A linear-phase tt/2 pulse can be re- 
focused with a gradient reversal, and a linear-phase 7T pulse can 
be used as a spin-echo pulse. The additional constraint of re- 
quiring linear phase has a cost in slice profile. If this constraint 
is eliminated better slice-profiles can be obtained. We present 
examples of minimum-phase and maximum-phase (reversed 
minimum-phase) pulses that have significantly better slice pro- 
files than the corresponding linear-phase pulses. These are use- 
ful when slice-profile phase is unimportant, such as inversion 
or saturation pulses. 

II. The Shinnar-Le Roux Algorithm 

The Bloch equation reduces to a rotation if relaxation effects 
are neglected. The solution to the Bloch equation is then the 
initial magnetization multiplied by a 3 X 3 orthonormal matrix. 
This rotation can also be represented by a 2 x 2 unitary matrix, 
the so-called spin-domain representation. The SLR algorithm is 
based on a discrete approximation to the spin domain version 
of the Bloch equation. This further simplifies the solution of the 
Bloch equation to the design of two polynomials. RF pulse de- 
sign becomes a polynomial design problem, which can be solved 
using well-known digital filter design algorithms. 

A. Spin Domain Bloch Equation 

The Bloch equation, neglecting relaxation, is given by 

k\ / o Gx -a,. A /a/A 

My ) =7 -Gx 0 If M v J 

M z / \ B Kv 0 / \M z J 

where fl, = B { x + iB Xy is the applied time varying RF field, 
and G is the gradient amplitude. The gradient is assumed to lie 
along the +.r axis. The solution to this equation is a rotation 

M(T) = RM(Q) 

where R is a 3 x 3 orthonormal matrix. If a sequence of n 
pulses has been applied, the solution is the product of the in- 
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dividual rotations /?, 



These rotations can also be represented by 2 x 2 unitary ma- 
trices 



■C "!*.) 



where a and 0 are the Cayley-KIein parameters 
a — cos 0/2 — in z sin 0/2 
j9 = -i(n, + in y ) sin 0/2. 



(1) 
(2) 



The vector n is the axis of the rotation and 0 is the rotation 
angle. These spin matrices multiply just as the 3 x 3 matrices 
do. The Cayley-KIein parameters satisfy the constraint 



(3) 



The spin matrices can be used to calculate the effect of an RF 
pulse in the presence of a constant gradient. We assume that the 
RF is piece-wise constant (illustrated in Fig. 1). This is a good 
approximation to the way that RF waveforms are generated on 
commercial imaging systems. If the jlh RF sample has a com- 
plex amplitude B t J and a duration At, the rotation parameters 
are 



0, = -yAty/\B x j\ 2 + (Gxf 
yAt 

nj = W\ 



(*,.,.,, a,. vj , Gx). 



The Cayley-KIein parameters for the jlh interval are then 
a s = cos 0 ; /2 - in z j sin 0,/2 
b j = -'K.y + in yJ ) sin 0,-/2. 

The spin-domain rotation matrix is then 



Qj = 



and the total rotation produced by the pulse is 
<2= 0.0,-1 ■•• C- 



(4) 



Given the spin-domain representation of the rotation pro- 
duced by a pulse, we can calculate the 3 x 3 rotation matrix by 
inverting the expressions for the Cayley-KIein parameters. 
Usually, though, the quantity of interest is the magnetization 
produced given some initial condition. These are calculated 
easily using the expression (adapted from fl7]) 



M*. \ 


/ («*) 2 




■ -(&*? 


m; J 








to \ 











where M xy = M x + iM y . 

There are a number of useful special cases. If the initial mag- 
netization is the equilibrium magnetization A/~ = (0, 0, M 0 ) T , 
the excitation slice profile is 



(5) 



B 
< 



1 
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1 

! 1 1 1 1 ill rfTTh-u 




if 



0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 I 

Time 

Fig. I. Piece-wise constant RF. This is a good approximation to the way 
RF waveforms arc generated on commercial imagers. 



and the inversion slice profile is 

M: = (aa* - W*)M 0 . (6) 

For the spin-echo profile the initial magnetization is M~ = 
(AC, A/~.*, 0) 7 . The spin-echo profile is then 

M+ = (cc*) 2 M; v - (?M-\ (7) 

If the initial magnetization is along the +y axis, this becomes 

(8) 

A case of considerable practical interest is of a spin-echo 
pulse surrounded by crusher gradients. The crusher gradients 
suppress the parasitic FID's produced by the slice profile tran- 
sition regions. A single crusher gradient produces a rotation 
about the z axis of 



Mxy = /((a*) + (j8) 2 )M 0 . 



on matrix for this rotati 
0 \ 



The spin-domain rotation matrix for this rotation is 



The rotation produced by the spin-echo pulse surrounded by 
crusher gradients is 



R ( RR C = 



K 0 a*e* 

The spin-echo pulse profile from (7) is 

AC = ( a *fe»*">M- - (0?M-*. 

The crusher gradient leaves a large number of cycles of phase 
twist across an imaged voxel. If we integrate over a voxel the 
first term integrates to zero, leaving 



AC = -0 2 M~*. 



(9) 



Again, if the initial magnetization was along the +y axis, the 
final magnetization produced by the crushed spin-echo pulse is 



(10) 



Once the Cayley-KIein parameters have been determined for 
a selective excitation pulse it is easy to calculate any profile of 
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interest. The spin domain description is an extremely conve- 
nient representation for selective excitation pulses. 

B. State-Space Description of Spin-Domain Rotations 

The result of an RF pulse applied in the presence of a con- 
stant gradient field can be calculated by multiplying 2x2 
unitary matrices, as in (4). This matrix product has considerable 
redundancy due to the symmetry of the spin-matrices, allowing 
an even simpler representation. The state of the rotation at a 
given point in the pulse sequence may be represented by a 2 x 
1 complex vector. The matrix-matrix products are replaced by 
matrix-vector products. 

If we write out the sequence of matrix products of (4) 




V ^ J 

Multiplication of the matrices up to j produces the rotation Q r 
From the structure of Q s it is apparent that two of the four en- 
tries are redundant. The same information can be represented 
by either column of the rotation matrix. We will choose the first 
column, {<xj pj) r . The effect of the pulse may then be calculated 
by propagating this vector 




This is a state-space description of the RF pulse. Since the ini- 
tial state is no rotation, the initial condition can be found by 
substituting a rotation angle <t> - 0 into (1) and (2). 

(:)-(:)• 

C. Hard-Pulse Approximation and the Forward SLR 
Transform 

The key step in the development of the SLR algorithm is the 
so-called hard-pulse approximation [18]. With this approxima- 
tion the state-space description can be reduced to two polyno- 
mials. This mapping of an RF pulse into two complex 
polynomials we call the forward SLR transform. 

During an excitation pulse the magnetization is rotating about 
the vector sum of the RF field and the local gradient field. The 
basic idea of the hard-pulse approximation is that if the angle 
is small the rotation can be modeled by two sequential rotations. 
The first is free precession under the effect of the local gradient 
field by an angle -yGx&t. The second is the rotation about the 
applied RF vector by an angle -7#,Ar. 

If we write Q } as the product of the two sequential spin ma- 
trices, we get 



where 

C; = cos( 7 |*,.,|A//2) 

S, o S i n ( 7 |*,jAl/2) 

The state-space recursion for the Cayley-Klein parameters then 

becomes 

CMr?)(: ;•)(::) 

Defining 

Aj = z'% 

Bj = z"% (12) 
the recursion may be reduced to 

(flj ~ {sj C y )(o z-')Gj-./ 

Given that A 0 = a 0 = 1 and B 0 - 0 O = °» tne fi rst tw0 states 

are 

CK:) 

/A 2 \ = /c 2 c, -s 2 *s,z-'\ 

\B Z ) \S 2 C X + C 2 S x z~ l I 

The Cayley-Klein parameters at the nth time step are (n - 1) 
order polynomials in z~ 1 

n - I 

A„(z) = S OjZ-' 

j=Q 

B a (z) = .2 bjz- J 

where 

z" 1 = e - hC '*'. 

The representation of the rotation produced by a selective ex- 
citation pulse has been reduced from a product of n 3 X 3 ma- 
trices to two (n - 1 )-order polynomials. 

The recursion given in (13) with the coefficients given in (1 1) 
maps an RF pulse into two complex polynomials A„(z) 
and B n (z). This mapping is the forward SLR transform. In the 
small-tip-angle case the SLR transform reduces to the familiar 

trans form. 

D. Inverse SLR Transform 

Given two related polynomials A„(z) and B„(z) y the SLR 
transform can be inverted to calculate the RF pulse that pro- 
duces these polynomials. This is called the inverse SLR trans- 
form. It is this inverse transform that reduces RF pulse design 
to polynomial design. 

The two (n - 1) order polynomials A„(z) and B„{z) must 
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satisfy the frequency response (slice profile) amplitude con- 
straint 

|^n(Of + \B„(Z)\ 2 = 1 (U) 

for all complex z such that | z \ = 1 . This is (3) with the substi- 
tutions indicated by (12). With this constraint the two polyno- 
mials A n (z) and B n (z) are valid representations of a rotation at 
any position x and z - e iyG * Al . 

The inverse SLR transform is found by the inversion of the 
state-space recursion. The forward recursion is 



(:)-(; t 



The inverse recursion is found by multiplying both sides by the 
inverse of the incremental rotation matrix, producing 



CjAj + Sfc 
z{SjAj + CjBj) 



(15) 



Since the incremental rotation and its inverse are both unitary 
matrices, the constraint of (14) is preserved, and the reduced- 
order polynomials are also valid representations of a rotation 
for all \z\ = 1. 

At any stage of the backward recursion we know the coeffi- 
cients of Aj(z) and Bj{z). Since Aj_ x {z) and Bj_ t (z) are lower 
order polynomials, the leading term in Aj_ } (z) and the low- 
order term in Bj. ,(z) must drop out 

qAjj^ =0 

-SjA JtQ + CjB jt0 = 0. 

These two equations are equivalent, as may be seen by expand- 
ing (14) as a polynomial and noting that all except the constant 
term are zero. The high-order term is 



-i^/o + B, 



With this relation either equation may be derived from the other. 

Here we choose the second equation to determine the RF. 
The ratio of the low-order terms is 



ie i9 > sin 0,/2 
cos <t>j/2 



where <fy is the tip angle produced by this segment of the pulse, 
and $j is the RF phase. The rotation angle is then 



<t>j = 2 tan 



and the RF phase is 

$j = *(-tf, i0 /4,„). 
The RF waveform is then 



(16) 



This expression for the RF coupled with the backward recursion 
of (15) constitute the inverse SLR transform. 

Hence, there is a unique invertible transform relationship be- 
tween an RF pulse and two polynomials A„{z) and B n (z) 

*S" {A n (z), B m (z)). 

This transform relationship makes RF pulse design equivalent 
to the design of two complex polynomials. The design of these 
polynomials is the subject of the rest of this paper. 

E. RF Pulse Design by Polynomial Design 

There are several different approaches that can be taken in 
designing the polynomials A n (z) and B„(z). One approach [1], 
[2 J, 16], (7J, [9] is to first approximate the magnetization com- 
ponent of interest with a polynomial, and then factor it accord- 
ing to the type of slice profile to solve for A n (z) and B„(z). 
The approach we will take here [7]— [9] uses the fact that 
B„(e nC **') is proportional to the sine of half the tip angle at 
position x (2). B n (z) is designed to optimally approximate the 
ideal slice profile. A„(z) is then calculated to be consistent with 
B„(z), subject to the additional constraint that the resulting RF 
pulse have minimum energy. Once A n (z) and B„(z) have been 
determined the RF pulse is found by the inverse SLR transform. 

For the sake of concreteness assume we are designing a x/2 
pulse about the jr-axis. In this case the ideal B,(z) polynomial 
evaluated at e' yCx *' is 

B,( e '-< c '") = J(«, + w.v) sin #(x)/2 

= i sin 0(jc)/2 

since n - ( 1, 0, 0) T . In-slice this is 

Bt{e h<UA ') = i sin x/4 = /V2/2 

while out-of-slice this is 

B f (e iyCxAt ) = /sinO = 0. 

Given this ideal B,{z) we wish to find a polynomial approxi- 
mation. There are many possible methods for designing this 
polynomial. One powerful method is the Parks-McClellan (PM) 
algorithm for the design of linear-phase finite impulse response 
(FIR) digital filters [19]. The algorithm requires the specifica- 
tion on the edges of the in-slice and out-of-slice regions, and 
the relative ripples in each. These can be calculated fairly ex- 
actly using the pulse design trade-off relationships that will be 
presented below. The ideal B,(z) and a polynomial approxi- 
mation both evaluated along the unit circle z = e iyCxA ' are plot- 
ted in Fig. 2. 

Given the B n (z) polynomial we need to find a consistent 
A n (z) polynomial. The constraint aa* + 00* = 1 means that 
the magnitude of the A n (z) polynomial is 



\A„{z)\ = Vl -B h {z)b:(z) (17) 

again evaluated along the unit circle z = e nGx *\ An additional 
constraint is required to make the choice of A n {z) unique. One 
possibility is to choose A„(z) to be a minimum-phase polyno- 
mial. As was mentioned in [7], a minimum phase A n {z) results 
in a minimum energy RF pulse. 

The minimum-phase A n (z) is easily found by noting that a 
minimum-phase polynomial is an analytic signal [20]. Analytic 
signals have the property that their log-magnitude and phase are 
a Hilbert transform pair. Hence, the minimum-phase A n (z) is 

- |-4„(z)|exp[i3C{log|4,(z)|}]. (18) 
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Fig. 2. 



Ideal Bt(z) and the polynomial approximation produced by the PM 
algorithm. 



A -0.8 -0.6 -0.4 -0.1 0 0.2 0.4 0.6 0.8 \ 

Position 

Fig. 3. Minimum-phase A„(z) corresponding to the B„(z) polynomial of 
Fig. 2. 



Here JC { • } is the Hilbert transform operator. The PM algo- 
rithm gives the polynomial coefficients of B n (z). B n (z) must 
then be evaluated on the unit circle, which can be done with a 
DFT. The DFT order should be significantly larger than n to 
reduce aliasing errors in subsequent calculations. Once B„(z) 
has been calculated M„(z)| can be found using (17). This is 
substituted into (18) to produce the minimum-phase polynomial 
A n (z) as evaluated along the unit circle. The polynomial coef- 
ficients are then found by an inverse DFT. Fig. 3 is a plot of 
the minimum-phase A n (z) as evaluated along the unit circle 
corresponding to the B n (z) of Fig. 2. 

The reason the minimum-phase A„(z) corresponds to a min- 
imum RF power is seen by examining the a 0 term of A n (z). 
From the forward recursion given in (13) it can be seen that a Q 
has the simple form 



fl 0 = C n Cn-\ 



C 2 C, 



the product of the half angle cosines for all the incremental RF 
rotations. These incremental rotations are all small. Using the 
small angle approximation cos 0/2 = 1 - 0 2 / 8 tne product can 
be written 

o 0 = (1 - #/8)(l - fli.,/8) • • • (1 - 0|/8)(1 - fl/8) 
where Bj = y\B lnj \At. After multiplication the first two terms 



a 0 = 1 - 1 S B 2 

= 1 -{(y*,f Z IB,/. 

j = 0 



(19) 



The second term is proportional to the pulse's average RF 
power. Hence, the pulse with minimum RF power has the max- 
imum a 0 . The polynomial with the maximum a 0 is the mini- 
mum-phase polynomial [20], [21]. 

The RF waveform computed by the inverse SLR transform 
of A„(z) and B n (z) is plotted in Fig. 4, and its slice profile in 
Fig. 5. Note that the ripples in-slice and out-of-slice are of con- 
stant amplitude, and that the slice is very well refocused. The 
M x ripple is larger in-slice because ten times more ripple was 
allowed when the B n (z) polynomial was designed. The reason 
for this, and for choosing the width and location of the transi- 
tion bands will be the subject of the next two sections. 



E 
< 




Fig. 4. SLR x/2 pulse corresponding to the 8„{z) and A„(z) polynomials 
of Figs. 2 and 3. 
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Fig. 5. Slice profile for the pulse plotted in Fig. 4. 

III. Pulse Design Parameter Relations 

The SLR algorithm reduces RF pulse design to the design of 
a single polynomial B n {z). The Parks-McClellan algorithm 
produces an optimal polynomial in the Chebyshev sense of min- 
imizing the maximum ripple. As inputs it requires that the stop- 
band and passband edges be specified, and that the relative 
ripple amplitudes in each of these bands be given. Pulse design 
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then consists of balancing the conflicting requirements for rip- 
ple, transition bandwidth, pulse length, slice width, and RF 
power. 

The efficiency of the SLR algorithm allows RF pulses to be 
designed iteratively. The calculation of a 128 point 7r/2 pulse 
*ukes about 10 s on a Sun 3/50. Slice profile parameters can be 
iteratively optimized until a satisfactory design is achieved. 

There is a better alternative. The characteristics of optimal 
FIR filters have been studied in detail [19). Empirical relation- 
ships have been derived for the various filter design parameters. 
With some modification these relationships can be used to pre- 
dict the attainable performance of RF pulses. Once a desired 
profile has been determined, the corresponding RF pulse is eas- 
ily and rapidly computed. 

The PM algorithm produces a linear-phase FIR filter. A lin- 
ear-phase ir/2 pulse can be refocused by a gradient reversal, 
and a linear-phase t pulse can be used as a spin-echo pulse. 
Optimal minimum-phase filters can be designed using the 
method described in [22], [23]. These have better performance, 
but at the cost of nonlinear phase response. When phase does 
not matter, as for selective inversion or saturation, the pulses 
based on minimum-phase filters will produce better slice pro- 
files. Maximum-phase filters are also of interest. These are sim- 
ply reversed minimum-phase filters. 

A. Recasting the FIR Filter Parameter Relations 

The parameter trade-offs for FIR filter design must be cast in 
a different form to be useful for the design of RF pulses. In 
particular the number of samples in the RF pulse is not ulti- 
mately significant, while the time-bandwidth product of the RF 
pulse is very important since it relates to the required RF power. 

Fig. 6 illustrates the key parameters of the performance of an 
FIR filter. The parameter 6, is the amplitude of the passband 
ripple, 6 2 is the amplitude of the stopband ripple, F p is the pass- 
band edge, and F s is the stopband edge. The passband ampli- 
tude is unity, and all frequencies are normalized to the Nyquist 
frequency . 

The empirical relationship between these parameters for an 
optimal FIR filter, such as is designed by the PM algorithm, is 
given by [19] 

6 2 ) = (N - \)AF +/(«,, d 2 )(AF)\ 

The parameter AF is the transition width F p - F,. The term 
0o»(6|, 6 2 ) is an empirically derived performance measure of 
the filter. The term containing the /(5,, 6 2 ) factor is negligible 
for almost any practical RF pulse, and will be dropped. The 
sample spacing is At, so the length of the pulse is T = (N - 
I) At. The half-amplitude width of the filter is approximately 
the average of the stopband and passband widths 



B = 



If, + 2F, 
2At 



If we define the fractional transition width as 

F p ~ F s 
F P + F s 

the equation may be simplified to 

«2) = TBW. 

This is the key design equation. Using it, the slice width, tran- 
sition width, pulse duration, and ripple amplitudes may all be 



(20) 




- 0.5 F 



Fig. 6. Digital filter parameters. 

traded off analytically. This expression also allows the pulse 
designer to exactly calculate the required inputs for the Park- 
McCIellan algorithm to produce an RF pulse with a specific 
performance. 

The function D w is an empirically derived performance mea- 
sure for optimal FIR filters. For linear phase FIR filters it is 
given by [19] 

D-A*i, «z) - + <hL, + ^ + [o«L? + a*/., + a 6 ] 

(21) 

where L x - log, 0 6,, ^ = log l0 b 2 , and the a { 's are the coeffi- 
cients 

a, = 5.309 x 10" 3 a 4 = -2.66 x 10" 3 

a 2 = 7.114 x 10 -2 a 5 = -5.941 x 10" 1 

a 3 = -4.761 x 10" 1 a 6 = -4.278 x 10"'. 

Fig. 7 is a plot of the linear-phase D w as a function of 6 2 for a 
range of values for 5,. 

Minimum-phase FIR filters are easily designed using the 
method described in [22], [23], If a minimum-phase filter of 
length N with ripple amplitudes 6, and 6 2 is desired, the PM 
algorithm is first used to design a linear-phase filter of length 
2N - 1 with ripple amplitudes 26, and 62/2. If the frequency 
response of this linear-phase filter is biased up so as to be al- 
ways positive, the result is the magnitude response squared of 
a minimum-phase filter. The minimum-phase filter is then com- 
puted by taking the square root to get the magnitude response, 
and then computing the corresponding minimum-phase poly- 
nomial as in (18). This filter will be of length A/, and have ripple 
amplitudes 5, and <5 2 . 

This minimum-phase filter design method performs the same 
operations as required in [6] to design B„{z) by factoring M z . 
It provides a direct way to compute one such factorization. In 
this case the method we use here of initially designing B„(z) 
produces one of the same pulses that is obtained by initially 
designing M, and factoring. 

The same basic design relation of (20) still holds for mini- 
mum-phase filters provided the correct expression for D M is 
used. £)„ m for optimal minimum-phase equal-ripple filters is 
given in terms of that for linear-phase filters by the expression 



/V M (6„ &i) = 1^/(25,^/2). 



(22) 



For most filters we will be concerned with (ripples of a few 
percent or less), the difference between D^ j and is on the 
order of 10-30%. This means that the minimum-phase filter can 
have a 10-30% narrower transition band, which is not a tre- 
mendous difference. However, for the same transition width the 
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Stopband Ripple, 6 2 

Fig. 7. / as a function of out-of-slice ripple. Each plot corresponds to 
a different value of the in-slice ripple, starting at 0.001 for the top plot and 
going to 1.0 for the bottom plot. Plots corresponding to 0.001 , 0.0 1 , 0. 1, 
and 1 are solid, 0.002, 0.02, and 0.2 arc dashed, and 0.005, 0.05, and 0.5 
are dot-dashed. 



minimum-phase filter passband ripple can be smaller by a factor 
of four or more. This is a significant difference. 

The design relation given in (20) describes the B n (z) poly- 
nomial. Unfortunately the relationship is generally nonlinear 
between the parameters of the B n (z) polynomial and the rele- 
vant feature of the slice profile. The nonlinearity is due to the 
geometry of the relationship between the B„(z) polynomial and 
the parameter of interest. Its effect is primarily in ripple ampli- 
tudes. We define the effective ripple amplitudes 5\ and 6 e 2 as the 
ripple actually produced in the magnetization component of in- 
terest relative to the maximum value of the profile. For exam- 
ple, for a tt/2 pulse 5' and 6 2 are the in-slice and out-of-slice 
ripples in M xy relative to M 0 . The 6, and 6 2 specified in the PM 
algorithm result in a slice profile with effective ripple ampli- 
tudes 6' and 6 2 . If the relationship between the specified poly- 
nomial ripple and effective slice profile ripple is known, we can 
calculate the ripple amplitudes to specify (6, and 6 2 ) in order 
to get the effective ripple amplitudes that we would like (6' and 
62). We will derive these relationships for five different types 
of pulses: small-tip-angle and tt/2 excitation pulses, inversion 
pulses, crushed spin-echo pulses, and saturation pulses. 

B. Small -Tip- Angle Pulses 

The small-tip-angle case is the simplest. Here the initial mag- 
netization lies along the z axis, and the parameter of interest is 
the transverse magnetization. In the small-tip-angle regime the 
RF pulse is proportional to the B n (z) polynomial, and the slice 
profile is proportional to the Fourier transform of the pulse. No 
compensation is required, so the specified and effective ripple 
amplitudes are the same 

6 2 = 6$. 

C. tt/2 Pulses 

The next case is of x/2 selective excitation pulses. Here there 
is considerable geometric distortion between the ripple as spec- 
ified by the B n (z) polynomial and the ripple in the resulting 






-±6 e 2 
(b) 




(d) 

Fig. 8. Ripple geometry for r/2 pulse. The out-of-slice ripple geometry 
is illustrated in (a) and (b). (a) is the half-angle or spin-domain geometry, 
while (b) is the full-angle or magnetization space geometry. The in-slice 
half-angle and full-angle ripple geometries are shown in (c) and (d). 



slice profile. To calculate the effective 6\ and P 2 we will assume 
that the B n {z) polynomial evaluated along the unit circle is 

B n (e iyCxA ') » sin 0(jc)/2 

where <p(x) is the tip angle at a position x. 

The out-of-slice ripple is in the small-tip-angle regime. From 
Fig. 8(a) it can be seen that for z corresponding to a ripple peak 

5 2 sin 7 a sin <p(jr)/2. 
4 

The sin( ir/4) factor scales the polynomial to the proper in-slice 
value for a x/2 pulse. The effective ripple is illustrated in Fig. 
8(b), and can be calculated as 

63 = sin <p(jc) 

= 2 sin <p(jc)/2 

V2 

= V262. 

The polynomial ripple that corresponds to a given effective rip- 
ple is then 

In-slice the C?„U) polynomial is bounded by (1 ± 6,) 
sin(ir/4). At an in-slice ripple peak 

(1 ± 6,) sin j = sin <p(*)/2. 

This is illustrated in Fig. 8(c). From the figure we deduce 
<t>(x) a tt/2 ± 26,. 
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The effective ripple is the variation in M xy produced by the 26, 
ripple in <f>(x). The effective ripple 6\, shown in Fig. 8(d), can 
then be approximated as 

5\ ^ 1 - cos 26, 

- 1 - (1 -i(26,) 2 ) 
= 26?. 

Note that the ripple is no longer symmetric about the desired 
value. The filter passband had a tolerance of ±6, while the slice 
profile has a tolerance of +0, -26?. The polynomial ripple cor- 
responding to a given effective ripple is 

6, = >/6t72. 

D. Inversion Pulses 

For inversion pulses the parameter of interest is M z : The ma- 
jor change in the analysis for v pulses is that the B„(z) poly- 
nomial is constrained to be less than one. The polynomial 
produced by the PM algorithm has a value of 1 ±6, in the 
passband. This is scaled so that the passband is bounded by ( 1 
- «.) ± 

From (6) the inversion pulse slice profile is given by 

M* = (aa* - 00*)M O . 

Using the constraint that aa* + 00* = 1 , this equation can be 
reduced to a function of 0 only 

Ml = (1 - 200*)M O . (23) 

An out-of-slice ripple of 6 2 in the B n (z) polynomial produces a 
ripple 

in the out-of-slice M z profile. The polynomial ripple corre- 
sponding to a given effective ripple is then 

5 2 = V6J72. 

In-slice B„(z) goes from (1 - 25, ) to 1 . The maximum in-slice 
ripple occurs when B n (z) is a minimum 

M, + « (1 - 2(1 - 26 1 ) 2 )A/ 0 
= (1 - 2 + 86, - S6])M 0 
= -I + 85,. 
The effective ripple is then 

Si = 86, 

and the polynomial ripple for a given effective ripple is 

6, = at/8. 

E. Crushed Spin-Echo Pulses 

The case of a spin-echo pulse surrounded by dephasing, or 
crusher, gradients is of considerable practical interest. The 
magnetization is initially in the transverse plane, and the param- 
eter of interest is the resulting M xy . 

From (10) the transverse magnetization produced by a crushed 
spin-echo RF pulse is 

M xy = i» 2 M 0 
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for an initial magnetization is along the + y axis. The effective 
out-of-slice ripple is then 

which may be inverted to give 

6 2 = V^. 

The B n {z) polynomial is scaled to be less than one, as it was 
for inversion pulses. In-slice the polynomial is bounded by ( 1 
- 6, ) ± 6,, so the effective ripple is 

« = 1 - (1 -26,) 2 
= 46,. 

Solving for 6, as a function of 6', 

6, = aj/4. 

F. Saturation Pulses 

Saturation technically means the destruction of the sample's 
net magnetization. It is then no longer able to produce a signal 
until the magnetization recovers. What we are actually referring 
to here are x/2 excitation pulses followed by a large dephasing 
gradient. This suppresses the signal produced by the magneti- 
zation in a slice, but does not destroy it. For example, it can 
still produce a stimulated echo. The term "suppression pulses'* 
would be more accurate, but the term "saturation pulses" is 
now well entrenched. 
The important magnetization component for saturation pulses 
* is M z . After the saturation pulse all transverse magnetization is 
assumed to be dephased. The saturated magnetization should be 
left with no M z so that subsequent excitation pulses will produce 
no signal. The unsaturated magnetization should be left with its 
full M z so that when subsequently excited it will produce its full 
signal. 

The saturation pulse profile is given by the same expression 
as for the inversion profile, given in (23) 

M; - (1 - 200*)Af o . 

The difference is in the scale factor for the 0 polynomial, and 
this changes the nature of the geometry. A saturation pulse is 
basically a 7r/2 pulse, so that the 0 polynomial is the PM poly- 
nomial scaled to sin tt/4 = V2/2. 

Out-of-slice the polynomial ripple is ±6 2 >/2/2, so that the 
out-of-slice M z ripple is 

Mt - 0 - 200* )M 0 

= (1 - 2(5 2 >/2/2)V 0 

= (1 - 5i)Af 0 . 

The effective ripple is then 

si -«t 

which can be inverted to give 

6 2 = V6J. 

i In-slice the polynomial ripple is bounded by ( 1 ± 6, ) V2/2, 
so the in-slice M z is bounded by 

I M: = (1 - 200*)M O 

- (1 " (1 ±5,) 2 )M 0 

= + 26,A/ 0 . 



O 
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TABLE I 

Parameter Relations for Various RF Pulses 



Case 5, 5 2 



Small-Tip 


5; 

V^72 


k 


w/2 


Inversion 


s;/8 




Spin-Echo 






Saturation 







The effective ripple is then 

«? = 25, 

so the polynomial ripple for a given effective ripple is 
«. = «T/2. 

G. Summary 

This section has been largely concerned with determining the 
relationship between the ripple in B n (z) and the effective ripple 
in the magnetic resonance parameter of interest. For conve- 
nience the key results are summarized in Table I. 

IV. Design Examples 

We next present several examples of how the results of the 
previous section are used to design pulses. The design relations 
first allow the pulse designer to trade off parameter values an- 
alytically. Then when an acceptable set of parameters has been 
selected, the relations tell exactly how to design the pulse that 
meets the predicted performance. 

A. Linear-Phase n /2 Pulses 

The first example is the design of slice-selective x/2 pulses. 
These are a principal component of most imaging pulse se- 
quences. This case will be examined in some detail both to il- 
lustrate the design method, and to show the nature of the design 
trade-offs. 

Usually we would like the excited slice to be refocused by a 
reversal of the slice select gradient. This can be assured by using 
a linear-phase B„(z) polynomial. Most of the pulses we will 
present in this section will be linear-phase, and hence refocus- 
able. Refocusing is not always required, though. An example 
is a field-of-view restriction pulse, like the slab-select pulse used 
in 3-D imaging. The phase profile of the slab is resolved by the 
imaging gradients. The phase is then suppressed by magnitude 
detection. When refocusing is not required significant advan- 
tages can be obtained by using a minimum-phase pulse. These 
include improved slice profile and reduced flow artifacts. An 
example of such a pulse is presented in the next section. 

We would like to produce a 0.5 cm thick slice with a pulse 
duration of 4.0 ms. If the gradient system can produce a max- 
imum gradient strength of 1 G/cm, the bandwidth of the pulse 
must be 

B = yGAx 

= (4.26 kHz/G)(l G/cm)(0.5 cm) 
= 2.13 kHz. 
The time-bandwidth product is then 

TB = (4 ms) (2.13 kHz) 
= 8.52. 



For convenience we will assume that the time-bandwidth prod- 
uct we are designing to is 8. The resulting pulses will produce 
0.47 cm slices at full gradient amplitude, or will produce 0.5 
cm slices at a gradient amplitude of 0.94 G/cm. 

The remaining design parameters are the in-slice and out-of- 
slice ripple amplitudes, and the width of the transition band. 
Choosing any two of these determines the third. Here we choose 
the two ripple amplitudes and then solve for the required frac- 
tional transition width using (20). We could just as well choose 
the transition width and use (20) to solve for and then pick 
two ripple amplitudes that will produce that D w . 

We also have to decide how many samples to use for the RF 
pulse. Fewer samples result in faster execution of the pulse al- 
gorithm. It also produces larger tip angles for each sample in 
the RF pulse, and this can result in the hard-pulse approxima- 
tion being violated. The number of points only has to be suffi- 
cient so that the slice-profile ripples are greater than the errors 
introduced by »Je hard-pulse approximation. Higher quality slice 
profiles require finer sampling. The effect of inadequate sam- 
pling is usually most noticable as nonequal out-of-slice ripple. 
For the examples presented in this section 65 samples are used 
for the RF waveform, which results in a maximum tip angle of 
12° per sample for pulses with a time-bandwidth product of 8. 

The other consideration for choosing the number of samples 
is that the sampling of the RF waveform will also excite side- 
lobes at multiples of the RF sampling frequency. The location 
and shape of these sidelobes can be calculated by a Fourier 
analysis since the sidelobe tip angles will be small. These side- 
lobes may be pushed further out either by interpolating the RF 
waveform as in [6], or by interpolating the B n (z) polynomial 
before the inverse SLR transform. The latter approach has the 
added benefit of reducing the errors in the hard-pulse approxi- 
mation, resulting in more accurate RF pulses. Increasing the 
filter order used for B„(z) is effectively just an interpolation 
method. 

Solving (20) for the fractional transition width 



As a first example let the in-slice ripple and the out-of-slice 
ripple both be 1 %. Using the expressions from Table I the re- 
quired polynomial ripples are 5, = 7.1% and 6 2 = 0.71%. 
These values are then substituted into the linear-phase 
expression (21) to give 

Dm j = 1.44. 
The fractional transition width is then 
1 44 

W= ±2 = 0.18. 

o 

The PM algorithm requires the location of the edges of the stop 
and passbands relative to the sampling frequency. If the RF 
pulse has 65 samples over 4 ms, the sampling frequency is 16 
kHz. The slice is 2 kHz in width, so it extends from ±1 kHz. 
The transition band is 

BW = (2kHz)(0.l8) - 0.36 kHz. 

This puts the passband edge at ( 1 - 0.36/2) = 0.82 kHz and 
the stopband edge at ( 1 + 0.36/2) =1.18 kHz. Normalizing 
by the sampling frequency 

F p = 0.82/16 = 0.05125 

F s = 1.18/16 = 0,07375. 



I. II 
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Fig. 9. The performance of the 1% ripple v/2 pulse. Plot (a) shows \M xy \ t 
while (b) is the same plot on a log scale. This shows that the out-of-slice 
ripple is close to the design value of 1 %. Plot (c) is of 1 - log, 0 | M ty \ 
which shows the in-slice ripple. It is also close to its 1 % design value. Plot 
(d) shows that the transition band is approximately the specified 360 Hz, 
and that the half amplitude point falls close to the intended 1 kHz. 



The only remaining input for the PM algorithm is the weighting 
of the in-slice and out-of-slice ripples. The in-slice ripple should 
be weighted by 1/6, and the out-of-slice ripple by 1 /5 2 - 

The RF pulse that results from these specifications is that 
plotted previously in Fig. 4. The performance of this RF pulse 
is shown in Fig. 9. Fig. 9(a) is the absolute value of the slice 
profile. Fig. 9(b) is the same plot on a log scale. The out-of- 
slice ripple is close to the design value of 1%. Fig. 9(c) is a 
plot of 1 - log, 0 \M xy \. This shows that the in-slice ripple is 
also close to its design value of 1%, although again it is ex- 
ceeded slightly. Fig. 9(d) shows that the transition width is ap- 
proximately the desired value of 360 Hz, and that the half 
amplitude point falls close to the intended 1 kHz. The transition 
width can be increased slightly if it is critical to achieve the 
specified ripple amplitudes exactly. 

Figs. 10 and 1 1 compare this SLR pulse to a Hamming win- 
dowed sine pulse, such as is commonly used on commercial 
imaging systems. Although the RF pulses are similar in char- 
acter there are several difFerences in the slice profiles. The SLR 
pulse has uniform ripple both in-slice and out-of-slice. The 
Hamming windowed sine has very little ripple in-slice, and only 
one significant sidelobe out-of-slice. The amplitude of this side- 
lobe is about 5%. The transition width of the windowed sine is 
50% wider than that of the SLR pulse. The SLR pulse would 
be preferable in cases where adjacent slices are being excited, 
such as multislice imaging. Both the absence of the sidelobe 
and the sharper profile reduce crosstalk between slices. 

In the previous example we simply chose the in-slice and out- 
of-slice ripple to be 1%. This determined the transition width. 
An important design trade-off is the relationship between the 
ripple amplitudes and the transition width. A large reduction in 
ripple amplitude can be obtained with a relatively small in- 
crease in transition width. Fig. 12 is a plot of two SLR pulses 
designed with 5% and 0.2% in-slice and out-of-slice ripple am- 
plitudes. Fig. 13 is a plot of the slice profiles produced by these 
pulses. The ripple amplitude has been reduced by a factor of 25 
at a cost of a factor of 2.5 in transition width. This overstates 
the difference in transition widths. If the same 5% ripple criteria 




Fig. 10. Comparison of the SLR 1% ripple pulse (solid line) to a Ham- 
ming windowed sine pulse (dashed line). 
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Fig. II. Slice profiles produced by the SLR 1 % ripple pulse (solid line) 
and the Hamming windowed sine pulse (dashed line). 
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Fig. 12. SLR pulses with ripple amplitudes of 5% (solid line) and 0.2% 
(dashed line). 

are used to measure both transition widths, the increase is less 
than a factor of 2. A very large reduction in ripple is obtained 
with a moderate reduction in slice sharpness. 

B. Minimum- Phase ir /2 Pulses 

Linear-phase pulses are not always required or desirable. One 
example is a field-of-view restriction pulse. Here the slice will 
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Fig. 13. Slice profiles produced by ihe 5% ripple (solid line) and 0.2% 
ripple (dashed line) SLR pulses. 
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Fig. 15. Refocused minimum-phase slice profile. The refocusing gradient 
lobe is 0.16 times the slice-select gradient area. 



0.2 




0 0.5 1 1.3 2 2.5 3~ 3.5 4 

Time, ms 

Fig. 14. Minimum-phase ir/2 pulse. 



be resolved by the imaging gradients, so the phase across the 
slice is unimportant provided it is small across a voxel. Nishi- 
mura et al. showed that an off-center RF pulse produces signif- 
icantly less flow dephasing than a linear-phase pulse, with only 
a minor penalty in slice profile [24]. A minimum-phase pulse 
obtains a similar reduction in flow dephasing along with an im- 
provement in slice profile. 

The minimum-phase pulse requires a maximum-phase filter 
because the filter phase is relative to its first sample, while the 
phase produced by the pulse is relative to its last sample. The 
maximum-phase filter is simply the minimum-phase filter re- 
versed. 

The minimum-phase pulse will have essentially the same in- 
tegrated RF power as the corresponding linear-phase pulse. Both 
have similar \ B n (z) \ profiles as evaluated along the unit circle, 
so both will have essentially the same A„(z) polynomials. From 
(19), both pulses will have the same integrated RF power. 

The higher performance of the minimum-phase filters can be 
used to improve any combination of the slice-profile parame- 
ters. If we again specify in-slice and out-of-slice ripple ampli- 
tudes at 1% the minimum-phase M may be calculated from 
(22) to be 1 .26. This is smaller than the 1 .44 for the linear- 
phase pulse. The reduction in D<» can be used to reduce any 
term in the product TBW. For example, the transition width W 
can be reduced by 15% to 300 Hz. This is 60 Hz less than for 
the linear-phase pulse. The RF waveform and the simulated slice 
profile are plotted in Figs. 14 and 15. The slice profile is shown 



refocused by a gradient area of only 0.16 times the slice select 
gradient area. In comparison a linear-phase pulse requires a re- 
focusing gradient area of 0.5 times the slice select gradient area. 

Rather than reducing one of the terms in the TBW product, 
we could also choose smaller ripple amplitudes so that Z>„. OT is 
the same as D mJ for the 1% ripple linear-phase pulse. Two 
possibilities are reducing the out-of-slice ripple by a factor of 
two, or reducing the in-slice ripple by a factor of five. 

C. 7T Inversion and Spin-Echo Pulses 

The same basic design method can be used to design inver- 
sion and spin-echo pulses. The only difference is which expres- 
sion from Table 1 is used to compute the polynomial ripple 
required to produce the desired effective ripple. The computed 
polynomial ripple is then used to compute D^. 

The example presented in this section has the same band- 
width and duration as the t/2 pulse of the previous section, 2 
kHz and 4 ms, respectively. The time-bandwidth product is 
again 8. The design is of an inversion pulse with 1 % ripple 
in-slice and out-of-slice. Since we do not care about phase for 
an inversion pulse we can use either a linear-phase or minimum- 
phase pulse. 

For the linear-phase case the corresponding D„ , is 2. The 
fractional transition width is 2/8 = 0.25, for a transition width 
of 500 Hz. The transition band then extends from 750 Hz to 
1.25 kHz. The number of points used in the design is 128 to 
ensure that the hard-pulse approximation is satisfied. The re- 
sulting pulse and its slice profile are shown as the solid lines in 
Figs. 16 and 17. The slice profile actually has slightly less rip- 
ple than specified. The in-slice and out-of-slice ripple are both 
about 0.8%. 

Note that the pulse has two large values at either end of the 
pulse. T^ese are due to the fact that a significant out-of-slice 
ripple in the spectrum of the pulse produces very little out-of- 
slice ripple in the M z component. This ripple allows a narrower 
transition width to be obtained. These are known in our group 
as "Conolly wings" because they appeared frequently in Con- 
olly's optimized inversion pulses [15]. 

For the minimum-phase case the D„ m is 1.36, which is sig- 
nificantly less than for the linear-phase pulse. The transition 
width is then reduced to 340 Hz. The minimum-phase pulse and 
its slice profile are plotted as the dashed lines in Figs. 16 and 
17. The ripple amplitudes are almost exactly the designed 1% 
values. 
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Fig. 17. Af. produced by the 1% ripple SLR inversion pulses of Fig. 16. 
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Fig. 18. Spin-echo profile produced by the linear-phase pulse of Fig. 16 
used without surrounding crusher gradients. 

The linear-phase pulse can also be used as a spin-echo pulse. 
The spin-echo slice profile produced by this pulse without sur- 
rounding crusher gradients is plotted in Fig. 18. The slice pro- 
file with surrounding crusher gradients is plotted in Fig. 19. In 
both cases the phase across the slice is very fiat. The ripple is 



less than for the inversion profile, as would be expected from a 
comparison of the entries in Table I. 

In [6] it is noted that an inversion pulse can have half the 
transition width of a spin-echo pulse of the same duration. Here 
this is reflected in the difference between D mJ for the refocusing 
pulse and m for the minimum-phase pulse. This difference 
depends on the number of passband zeros in the B„ ( z ) poly- 
nomial [22 J. For a very narrow slice with no passband zeros (a 
single lobe RF pulse) the minimum-phase and linear-phase 
pulses are identical. For a wide passband with many passband 
zeros the 2 to 1 ratio noted in [6] can be approached. For the 
examples shown in this section the minimum-phase inversion 
pulse has a transition width that is 340/500 = 0.68 times that 
of the linear-phase pulse. 

The design of x pulses is probably the most important appli- 
cation of the SLR algorithm. Fourier methods fail badly in the 
design of pulses with tip angles of this size. Optimization meth- 
ods are capable of designing good v pulses, but require consid- 
erably more time and computation. The SLR algorithm makes 
the design of tt pulses quick, and it analytically predicts the 
attainable slice profile. 

D. Saturation Pulses 

An ideal saturation pulse would leave no M z magnetization 
for the saturated spins. In practice the effectiveness of a satu- 
ration pulse is determined primarily by the RF uniformity across 
the subject. On our 1.5T GE Signa System the subject induced 
variation is on the order of 10%. The saturation pulse in-slice 
ripple should be somewhat less than this so that RF uniformity 
is the performance limiting factor. The out-of-slice ripple can 
easily be made quite small due to the square dependence of the 
M z ripple on the polynomial ripple. 

As an example we again consider a 4 ms, 2 kHz bandwidth 
pulse. Since the slice will not be refocused we do not need a 
linear-phase pulse. The excited magnetization will be sup- 
pressed by a subsequent dephasing gradient anyway, so it is 
advantageous for the excitation pulse to leave the slice with as 
much phase as possible. This is accomplished by using a max- 
imum-phase pulse, which is based on a minimum-phase B„{z). 

If we specify the in-slice ripple as 1% and the out-of-slice 
ripple as 0. 1 % , we calculate D„ to be about 1.36. This gives a 
fractional transition width of 0.17, or a transition width of 340 
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Fig. 20. Maximum-phase saturation pulse. 




Frequency, kHz 

Fig. 21. Log-magnitude slice profile for the saturation pulse of Fig. 20. 



Hz. The resulting pulse is shown in Fig. 20, and a log plot of 
the slice profile is shown in Fig. 21. 

o 

V. Conclusion 

This paper has outlined the Shinnar-Le Roux algorithm for 
the design of RF pulses. Beyond its inherent speed, perhaps its 
most useful feature is that slice profiles can be predicted ana- 
lytically. Slice profile parameters can be traded off before ever 
designing a pulse. This allows RF pulse design to be fully in- 
tegrated into pulse sequence design. A number of pulse designs 
were presented. These included the usual ir/2 and tt pulses used 
in most MR scanners, as well as a saturation pulse, and mini- 
mum-phase and maximum-phase pulses. 
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